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Abstract 



Lie symmetry group method is applied to study the boundary-layer equations for two-dimensional steady flow of an incom- 
|r^ ' , pressible, viscous fluid near a stagnation point at a heated stretching sheet placed in a porous medium equation. The symmetry 
.^ ■ group and its optimal system are given, and group invariant solutions associated to the symmetries are obtained. Finally the 
^s ' structure of the Lie algebra symmetries is determined. 
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1 Introduction 



> 

In fluid dynamics, a stagnation point is a point in a flow field where the local velocity of the fluid is zero. Stagnation 
C^ ' points exist at the surface of objects in the flow field, where the fluid is brought to rest by the object. The Bernoulli 
I — | , equation shows that the static pressure is highest when the velocity is zero and hence static pressure is at its 
f ^ ■ maximum value at stagnation points. This static pressure is called the stagnation pressure. The Bernoulli equation 
f~) applicable to incompressible flow shows that the stagnation pressure is equal to the dynamic pressure plus static 
t-H . pressure. Total pressure is also equal to dynamic pressure plus static pressure so, in incompressible flows, stagnation 
: pressure is equal to total pressure. (In compressible flows, stagnation pressure is also equal to total pressure providing 
. £_< ' the fluid entering the stagnation point is brought to rest is entropically.) 

x : 

?-h ' In physics and fluid mechanics, a boundary layer is that layer of fluid in the immediate vicinity of a bounding surface. 



In the Earth's atmosphere, the planetary boundary layer is the air layer near the ground affected by diurnal heat, 
moisture or momentum transfer to or from the surface. On an aircraft wing the boundary layer is the part of the 
flow close to the wing. The boundary layer effect occurs at the field region in which all changes occur in the flow 
pattern. The boundary layer distorts surrounding non- viscous flow. It is a phenomenon of viscous forces. This effect 
is related to the Reynolds number (In fluid mechanics and heat transfer, the Reynold's number is a dimensionless 
number that gives a measure of the ratio of inertial forces to viscous and, consequently, it quantifies the relative 
importance of these two types of forces for given flow conditions) . Laminar boundary layers come in various forms 
and can be loosely classified according to their structure and the circumstances under which they are created. The 
thin shear layer which develops on an oscillating body is an example of a Stokes boundary layer, whilst the Blasius 
boundary layer refers to the well-known similarity solution for the steady boundary layer attached to a flat plate 
held in an oncoming unidirectional flow. When a fluid rotates, viscous forces may be balanced by the Coriolis effect, 
rather than convective inertia, leading to the formation of an Ekman layer. Thermal boundary layers also exist in 
heat transfer. Multiple types of boundary layers can coexist near a surface simultaneously. The deduction of the 
boundary layer equations was perhaps one of the most important advances in fluid dynamics. Using an order of 
magnitude analysis, the well-known governing NavierStokes equations of viscous fluid flow can be greatly simplified 
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within the boundary layer. Notably, the characteristic of the partial differential equations (PDE) becomes parabolic, 
rather than the elliptical form of the full NavierStokes equations. This greatly simplifies the solution of the equations. 
By making the boundary layer approximation, the flow is divided into an inviscid portion (which is easy to solve by 
a number of methods) and the boundary layer, which is governed by an easier to solve PDE. 

Flow and heat transfer of an incompressible viscous fluid over a stretching sheet appear in several manufacturing 
processes of industry such as the extrusion of polymers, the cooling of metallic plates, the aerodynamic extrusion 
of plastic sheets, etc. In the glass industry, blowing, floating or spinning of fibres are processes, which involve the 
flow due to a stretching surface. Mahapatra and Gupta studied the steady two-dimensional stagnation-point flow 
of an incompressible viscous fluid over a flat deformable sheet when the sheet is stretched in its own plane with a 
velocity proportional to the distance from the stagnation-point. They concluded that, for a fluid of small kinematic 
viscosity, a boundary layer is formed when the stretching velocity is less than the free stream velocity and an inverted 
boundary layer is formed when the stretching velocity exceeds the free stream velocity. Temperature distribution in 
the boundary layer is determined when the surface is held at constant temperature giving the so called surface heat 
flux. In their analysis, they used the finite-differences scheme along with the Thomas algorithm to solve the resulting 
system of ordinary differential equations. 

This paper is concerned with the solution of steady two-dimensional stagnation point flow of an incompressible 
viscous fluid over a stretching sheet which is placed in a fluid saturated porous medium. Lie-group theory is applied 
to the equations of motion for determining symmetry reductions of partial differential equations. 

2 Lie Symmetries of the Equations 

A PDE with p— independent and q— dependent variables has a Lie point transformations 
Xi = xi + s£i(x, u) + 0(e 2 ), u a = u a + E(p a (x, u) + 0(e 2 ), 

for a = 1, ..., q. The action of the Lie group can be considered 
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v = ^&(x,u) 7 r- + 5> Q (z,u) 7 r- (1) 



on the total space of PDE (the space containing independent and dependent variables) . Furthermore, the character- 
istic of the vector field (I) is given by 



P a a 

Q Q (x,u (1) ) =(p a (x,u) -^2^i(x,u)-^-, 
and its n~ th prolongation is determined by 

i=\ % a=l JJ=j=0 J 

where <p J a — DjQ a + Y^d=i £i u ji- (Dj is the total derivative operator describes in (4)). 

The deduction of the boundary layer equations was perhaps one of the most important advances in fluid dynamics. 
Using an order of magnitude analysis, the well-known governing NavierStokes equations of viscous fluid flow can be 
greatly simplified within the boundary layer. Notably, the characteristic of the partial differential equations becomes 
parabolic, rather than the elliptical form of the full NavierStokes equations. This greatly simplifies the solution of 
the equations. By making the boundary layer approximation, the flow is divided into an inviscid portion (which is 
easy to solve by a number of methods) and the boundary layer, which is governed by an easier to solve PDE. The 



continuity and NavierStokes equations for a two-dimensional steady incompressible flow in Cartesian coordinates 
are given by 
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where U,V,P and T are smooth functions of (x,y), and also U and V are the velocity components along x-axis 
and y-axis respectively, v is the kinematic viscosity, k is the permeability of the porous medium, T is the fluid 
temperature, a is the coefficient of thermal diffusivity. The aim is to analysis the Lie point symmetry structure of 
the system of steady two-dimensional boundary-layer stagnation-point flow equations. 

Let us consider a one-parameter Lie group of infinitesimal transformations [x, y, U, V, P, T) given by 



x = x + e^(x, y, U, V, P, T) + 0(e 2 ), 
U = U + e m (x,y,U,V,P,T) + 0(e 2 ), 
P = P + e m {x,y 1 U 1 V 1 P 1 T) + 0(e 2 ) 1 



y = V + e&(x,y,U,V,P,T) + 0(e 2 ), 
V = V + e m (x,y, U, V,P,T) + 0(e 2 ), 
f = T + e m (x, y, U, V, P, T) + G(e 2 ), 



where e is the group parameter. Then one requires that this transformations leaves invariant the set of solutions of the 
system (2). This yields to the linear system of equations for the infinitesimals £i(x, y, U, V, P, T), £2(2;, y, U, V, P, T), 
r}i(x,y,U, V, P, T) , 7/2 (x, y, U, V, P,T),r]s(x,y, U, V, P, T) and 7/4 (x, y, U, V, P, T) . The Lie algebra of infinitesimal sym- 
metries is the set of vector fields in the form of 
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Iff = D x ((fl - ClU x - &Uy) + £lU xx + &U x y, 

d I(t°i - £,iU x - &U y ) + t;iU xxx + &U XX , 
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D x Dy(tf2 - i\V x - frVy) + (,\V xx y + &V xx y, 

D x (<p 3 ~ £lP x - &Py) + &P XX + &Pxy, 
Ar(^3 ~ £,\Px — &Py) + £,lPxxx + ^Pxxy, 
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Table 1 

Commutation relations of g 

[ , ] Vi v 2 v 3 v 4 

vi 

v 2 

v 3 v 3 

V 4 -V3 



where the operators D x and D y denote the total derivative with respect to x and t: 

D * = i +u °w +v °w +p *^P +T ^-> D » = ^ y +U »w +V «w +P ^ +T «wf---- (4) 

Using the invariance condition, i.e., applying the second prolongation v*- 2 ', (3), to system (2), and by solving the 
linear system 
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the following system of 112 determining equations yields: 
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The solution of the above system gives the following coefficients of the vector field v: 

£i = Ci, £ 2 = C 2 , r?i =0, ry 2 = 0, 773 = T, 774 = 0, 

where Ci and C2 are arbitrary constants, thus the Lie algebra $3 of the steady two-dimensional boundary-layer 
stagnation-point flow equations is spanned by the four vector fields 

d d d d 

which vi,v 2 and V3 are translation on x,u and T, V4 is scaling on T. The commutation relations between these 
vector fields is given by the table 1, where entry in row i and column j representing [vj,Vj]. 

The one-parameter groups Gj generated by the base of q are given in the following table. 

G!:(x + e,y,U,V,P,T), G 2 : (x,y + e,U,V,P,T), G 3 : (x,y,U,V,P,T + e), G 4 : (x,y, U, V,P,Te £ ). 
Since each group Gj is a symmetry group and if U = f(x, y), V — g(x, y),P = h(x, y) and T — r(x, y) are solutions 



of the system (2), so are the functions 

Ui= f(x + e,y), Ui = f(x,y + e), U 1 =f(x,y), U Y = f(x, y), 

Vi=g(x + e,y), W = g{x,y + e), V 1 =g(x,y), V 1 = g{x,y), 

Pi = h(x + e,y), Pi = h(x,y + s), Pi=h(x,y) 1 P 1 = h(x,y), 

T 1 = r{x + s,y), T 1 =r(x,y + e), T 1 =r{x,y)+e, T x = e~ s r(x,y), 

where e is a real number. Here we can find the general group of the symmetries by considering a general linear 
combination civi + • • • + C1V4 of the given vector fields. In particular if g is the action of the symmetry group near 
the identity, it can be represented in the form g = exp(e4V4) o • • • o exp(eivi). 

3 Symmetry reduction for steady two-dimensional boundary-layer stagnation-point flow equations 

Lie-group method is applicable to both linear and non-linear partial differential equations, which leads to similarity 
variables that may be used to reduce the number of independent variables in partial differential equations. By 
determining the transformation group under which a given partial differential equation is invariant, we can obtain 
information about the invariants and symmetries of that equation. 

The first advantage of symmetry group method is to construct new solutions from known solutions. Neither the first 
advantage nor the second will be investigated here, but symmetry group method will be applied to the Eq. (2) to 
be connected directly to some order differential equations. To do this, a particular linear combinations of infinitesi- 
mals are considered and their corresponding invariants are determined. The steady two-dimensional boundary-layer 
stagnation-point flow equations expressed in the coordinates (x,y), so to reduce this equation is to search for its 
form in specific coordinates. Those coordinates will be constructed by searching for independent invariants (r, s) 
corresponding to an infinitesimal generator. So using the chain rule, the expression of the equation in the new coor- 
dinate allows us to the reduced equation. Here we will obtain some invariant solutions with respect to symmetries. 
First we obtain the similarity variables for each term of the Lie algebra rj, then we use this method to reduced the 
PDE and find the invariant solutions. Here our system has two-independent and four-dependent variables, thus, the 
similarity transformations and invariant functions with invariant solutions are coming in the following table. 

vector field invariant function invariant transformations similarity transformations 

x = r,y = s, 
x = r,y = s 

translating of all function T is the similarity 

transformation 

scaling of all function T is the similarity 

transformation 

In the above table, the invariant solution with respect to vi and V2 obtained from translation on both dependent 
variables, and for V3 and V4 all invariant solutions are obtained by translating and scaling on function T. 

4 Optimal system of steady two-dimensional boundary-layer stagnation-point flow equations 

Let a system of differential equation A admitting the symmetry Lie group G,be given. Now G operates on the set of 
solutions S of A. Let s ■ G be the orbit of s, and H be a subgroup of G. Invariant H— solutions s £ S are characterized 
by equality s ■ S — {s}. If h € G is a transformation and s e S^then h- (s ■ H) = (h- s) ■ (hHh^ 1 ). Consequentlyevery 
invariant H— solution s transforms into an invariant hHhr 1 — solution (Proposition 3.6 of [5]). 

Therefore, different invariant solutions are found from similar subgroups of G. Thus, classification of invariant 
H— solutions is reduced to the problem of classification of subgroups of G,up to similarity. An optimal system of 
s— dimensional subgroups of G is a list of conjugacy inequivalent s— dimensional subgroups of G with the property 
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that any other subgroup is conjugate to precisely one subgroup in the list. Similarly, a list of s— dimensional subal- 
gebras forms an optimal system if every s— dimensional subalgebra of g is equivalent to a unique member of the list 
under some element of the adjoint representation: f) = Ad(g) • f). 

Let H and H be connected, s— dimensional Lie subgroups of the Lie group G with corresponding Lie subalgebras f) 
and t) of the Lie algebra g of G. Then H — gHg^ 1 are conjugate subgroups if and only f) = Ad(g) • f) are conjugate 
subalgebras (Proposition 3.7 of [5]). Thus, the problem of finding an optimal system of subgroups is equivalent to 
that of finding an optimal system of subalgebras, and so we concentrate on it. 

4-1 One- dimensional optimal system 

For one-dimensional subalgebras, the classification problem is essentially the same as the problem of classifying the 
orbits of the adjoint representation, since each one-dimensional subalgebra is determined by a nonzero vector in Lie 
algebra symmetries of steady two-dimensional boundary-layer stagnation-point flow equations and so to " simplify" 
it as much as possible. 

The adjoint action is given by the Lie series 



Ad(exp(evi)v,-) = v,- - e[vi, v,-] + -[v„ [vj, v,-]] - 



(5) 



where [v,, Vj] is the commutator for the Lie algebra, e is a parameter, and i,j — 1, • • ■ ,4. Let Ff : g — > g defined 
by vh Ad(exp(evi)v) is a linear map, for i = 1, • • • , 4. The matrices Mf of Ff, i = 1, • • ■ ,4, with respect to basis 
{vi,--- ,v 4 } are 
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by acting these matrices on a vector field v alternatively we can show that a one-dimensional optimal system of g is 
given by 



Xi = aiVi + a 2 V2, X 2 = aiVi + a 2 v 2 + a 3 v 3 , 
where a,-'s are real constants. 



(7) 



4-2 Two-dimensional optimal system 

Next step is to construct two-dimensional optimal system, i.e., classification of two-dimensional subalgebras of g. 
The process is by selecting one of the vectors in (7), say, any vector of (7). Let us consider Xi (or X 2 ). Corresponding 
to it, a vector field X = aiVi + • • • + 04 V4, where a^'s are smooth functions of (x, y, U, V, P, T) is chosen, so we must 
have 



[X 1 ,X] = \Xi+nX, 
the equation (8) leads us to the system 



(8) 



C l j k ajdk = A<2i + \xa.i 



(i= 1,2,3,4). 



(9) 



The solutions of the system (9), give one of the two-dimensional generator and the second generator is X\ or, X 2 
if selected. After the construction of all two-dimensional subalgebras, for every vector fields of (7), they need to be 
simplified by the action of (6) in the manner analogous to the way of one-dimensional optimal system. 



Table 2 

Commutation relations of g 



dimension 



{Q1V1+Q2V2} {/3lVl +/?2V2,/?3Vl +/?4V2 + V3,) (vi,V2,V3> (vi , V2, V3, V4} 

subalgebras (aivi + «2V2 + a 3 v3) (/Sivi + /32V2, V4) (vi,V2,V4> 

(V3,V 4 ) (V1,V 3 ,V3> 



(V2,V3, V4) 



Consequently the two-dimensional optimal system of g has three classes of g's members combinations such as 

1) /3lVi+/?2V2,/3 3 Vi+/?4V2 + V3, /?l/3 3 ^ foPi, (10) 

2) / 3iv 1 +/3 2 v 2 ,v 4 , ft 2 +ft 2 /0, 

3) v 3 ,v 4 . 

4--S Three-dimensional optimal system 

This system can be developed by the method of expansion of two-dimensional optimal system. For this take any 
two-dimensional subalgebras of (10), let us consider the first two vector fields of (10), and rename them Y\ and Y 2 , 
thus, we have a subalgebra with basis {yi, Y" 2 }, find a vector field Y = aivi + • • • + a 4 V4, where dj's are smooth 
functions of (x, y, U, V, P, T), such the triple {Yi, Y2, Y} generates a basis of a three-dimensional algebra. For that it 
is necessary an sufficient that the vector field Y satisfies the equations 

\Yi,Y} = XiY + fi 1 Y 1 + v 1 Y 2 , 

[Y 2 ,Y] = \ 2 Y + m Y 1 + v 2 Y 2 , (11) 

and following from (11), we obtain the system 

C l jk f3 J r a k = \ia i + iJ,il3 l r + is 1 pi, r = 1,2, s = 3,4, a = 1,2,3,4, (12) 

C i jk Pia k = \ 2 a i +H2Pi + V2PI, r=l,2, s = 3,4, a= 1,2,3,4. 

The solutions of system (12) is linearly independent of {Y\, Y 2 } and give a three-dimensional subalgebra. This process 
is used for the second two couple vector fields of (12). 

Consequently the three-dimensional optimal system of g is given by 

1) Vi, v 2 , v 3 , 2) vi, v 2 , v 4 , 

3) Vi, v 3 , v 4 , 4) v 2 , v 3 , v 4 . 

All previous calculations lead to the table (2) for the optimal system of q. 

5 Lie Algebra Structure 

q has a no any non-trivial Levi decomposition in the form of Q = t x jji, because q has no any non-trivial radical, 
i.e., if r be the radical of g, then g = r. 

If we want to integration an involuting distribution, the process decomposes into two steps: 

• integration of the evolutive distribution with symmetry Lie algebra g/t, and 

• integration on integral manifolds with symmetry algebra r. 



First, applying this procedure to the radical r we decompose the integration problem into two parts: the integration 
of the distribution with semisimple algebra q/x, then the integration of the restriction of distribution to the integral 
manifold with the solvable symmetry algebra r. 

The last step can be performed by quadratures. Moreover, every semisimple Lie algebra q/x is a direct sum of simple 
ones which are ideal in q/x. Thus, the Lie-Bianchi theorem reduces the integration problem ti evolutive distributions 
equipped with simple algebras of symmetries. Thus, integrating of system (2), become so much easy. 

Both q is a solvable, because if g^ — (vj,[vj,Vj]) = [q,q], we have gW = [q,q] = (vi,--- ,v 4 ), and q^ = 
[q^\q^] = ( v 3)j so, we have a chain of ideals gW D g*- 2 -* Z> {0}. According to the table of the commutators we can 
decompose g in the direct sum of g = R 2 © a(l), where a(l) is the Lie algebra of Afhne transformations group A(l). 

6 Conclusion 

In this article group classification of steady two-dimensional boundary-layer stagnation-point flow equations and the 
algebraic structure of the symmetry group is considered. Classification of r-dimensional subalgebra is determined 
by constructing r-dimensional optimal system. Some invariant objects are fined and the Lie algebra structure of 
symmetries is found. 
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